Field induced quantum-Hall ferromagnetism in suspended bilayer graphene 
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We have measured the magneto-resistance of freely suspended high-mobility bilayer graphene. 
For magnetic fields B > 1 T we observe the opening of a field induced gap at the charge neutrality 
point characterized by a diverging resistance. For higher fields the eight-fold degenerated lowest 
Landau level lifts completely. Both the sequence of this symmetry breaking and the strong 
transition of the gap-size point to a ferromagnetic nature of the insulating phase developing at the 
charge neutrality point. 



PACS numbers: 73.22.Pr,73.43.-f,71.70.Di,73.43.Qt 

I. INTRODUCTION 

The unique electronic properties of monolayer and bi- 
layer graphene makes them promising candidates for fu- 
ture applications in nanotechnology. Though (bilayer) 
graphene on a Si02-substrate can show a mobility up 
to 2 m 2 /VsPmuch cleaner and higher mobility samples 
are necessary in order to investigate its intrinsic proper- 
ties, and, in particular, electron interaction effects. Mo- 
bilities exceeding 10 m 2 /Vs can be obtained b y re mov- 
ing the SiC>2 substrate underneath the graphene^ or by 
depositing graphene on a boron nitride crystal.^ These 
high-mobility samples display new interaction-induced 
phenomena such as a fractional quantum Hall effect PS 
broken-symmetry states,^ a magnetic-field induced insu- 
lating phase,- and quantized conductance at zero mag- 
netic fieldP 

In the two-dimensional electron system of bilayer 
graphene (BLG) the application of a perpendicular mag- 
netic field results into an unconventional integer quan- 
tum Hall effect with plateaus at filling factors v — 
±4, ±8, ±12, ..PThe lowest Landau level is eight-fold 
degenerate, owing to spin, valley and layer-index de- 
grees of freedom. In standard BLG samples deposited 
on Si02, magnetic fields around 10 T are required to ob- 
serve fully quantized plateaus and the eight-fold degener- 
acy of the lowest Landau level is only lifted for the high- 
est quality samples at magnetic fields exceeding 20 TP^ 
At T the density of states in BLG does not vanish 
at the charge neutrality point, in contrast to single layer 
graphene, therefore, even arbitrarily weak interaction be- 
tween charge from conduction and valence band states 
will trigger excitonic instabilities which causes a variety 
of gapped statesPHUl 

In this paper, we present two-terminal magnetotrans- 
port experiments in suspended BLG at temperatures 
ranging from 1.3 K to 4.2 K and magnetic fields up to 
30 T. We observe a sudden gap opening at the CNP 
already for B > 1 T and the appearance of broken- 
symmetry states at filling factors v = ±1,±2,±3 for 



higher fields. Detailed investigation of the energy gap at 
filling factor v = reveals an exchange-interaction driven 
linear scaling at low magnetic fields, in agreement with 
earlier reported resultsP At high fields we observe the 
cross-over to a much smaller gap. This high field transi- 
tion and the appearance of broken symmetry states at v 
= 1,2,3 are consistent with t he formation of a quantum 
Hall ferromagnetic state! 1 ^ 1 ^ 



II. EXPERIMENTAL BACKGROUND 

We have prepared a suspended BLG sample using 
an acid free method.^ Following standard techniques 
we first exfoliated flakes from highly oriented pyrolytic 
graphite (HOPG) and deposited them on a Si/Si02 sub- 
strate covered with a 1.15 /im thick LOR- A resist layer. 
Bilayer flakes were then identified by their optical con- 
trast Subsequently, two electron beam lithography 
steps were performed in order to contact the flakes with 
Ti-Au contacts and to remove part of the LOR-A below 
the graphene flakes. The resulting device is freely sus- 
pended across a trench formed in the LOR-A with two 
metallic contacts on each side, see inset of Fig. [T] 

Carriers in the BLG sheet can be induced by applying 
a back-gate voltage Vg on the highly n-doped Si wafer. 
The geometrical gate capacitance is given by a combi- 
nation of the vacuum gap (1.15 /im) and Si02 substrate 
(0.5 /im). Using a serial capacitor model we calculate a 
gate capacitance of 7.2 aF//im~ 2 which directly relates 
the carrier concentration to Vg as n = a(Vct — Vcnp) 
with leverage factor a = 0.5 x 10 14 m _2 V~ 1 and a finite 
voltage of the the CNP of Vcnp = 1-2 V. In high mag- 
netic fields, the geometric capacitance increases due to 
the formation of edge states^ and a becomes dependent 
on B. Therefore, the exact values of capacitance were 
determined experimentally by identifying the filling fac- 
tors of quantized Hall plateaus in magnetic field, details 
can be found in the appendix. 

After mounting the devices were slowly cooled down to 
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4.2 K and current annealed^ by applying a DC bias cur- 
rent up to 3 mA. This local annealing resulted into the 
high quality sample with mobility /i « 10 m 2 /Vs at a 
charge carrier density n = 2xl0 11 cm~ 2 . The value of 
the mobility is calculated based on the dimension of sus- 
pended graphene before current annealing: 0.3 /im wide 
and 2.1 /im long. However, in the membrane the distri- 
bution of the temperature while current annealing is non 
homogenous,^ which most probably leads to the middle 
part of the membrane being annealed and non annealed 
regions close by the contacts. In this case the estimation 
of the mobility value based on geometrical dimensions 
might be not precise. We can also estimate the quality of 
obtained sample from the value of magnetic field at which 
the system enters the quantum Hall regime (B > 0.5 T). 
Assuming pB ^> 1 for QHE to existP, the observation 
applies a lower bound for the mobility of 2 m 2 /Vs. 
Measurements were performed with standard low- 
frequency lock-in techniques in two-probe geometry with 
an excitation current of 2 nA. 



III. RESULTS 

In Fig. [I] we show the data for the two-point resistance 
R of our suspended BLG device at B = T and B = 1 T 
as a function of Vq (top x-axis) and n (bottom x-axis), 
respectively. The two-probe resistance R is characterized 
by a magnetoresistance p xx — L/w ■ R xx with superim- 
posed Hall-resistance p xy , R — (L/w) ■ p xx + p xy . Here 
L/w « 6.7 is the aspect ratio of the device. The traces 
are corrected by phenomenological contact resistances 
(1 kfi on the electron-side and 1.7 k£! on hole-side) which 
were determined from a finite resistance background ob- 
served at high carrier concentrations; this background 
resistance increases by about a factor 2 in the range 
B = 0...30 T. These contact resistances most probably 
originate from in-series con nected non-annealed parts of 
the samplepS contact doping^SESland the finite resistance 
of the current leads. The sharp maximum at the CNP of 
the zero-field data already indicates the high electronic 
quality of the sample. At 1 T the resistance already ex- 
hibits fully quantized plateaus at filling factors v = 4 and 
a developing quantization at v = 8 and v = 12. The for- 
mation of these plateaus is caused by a quantization of 
Pxy = h/ve 2 and the associated zero minima in p xx when 
the Fermi energy lies between two Landau levels^ and 
confirms the high electronic mobility (p 3> 1/-B) of our 
device required to observe this unconventional quantum 
Hall effect. 

Additionally, as soon as a finite magnetic field is ap- 
plied, the resistance at the CNP, Rcnp, starts to di- 
verge. Whereas at zero magnetic field Rcnp is only very 
weakly temperature dependent and comparable to the 
resistance quantum, already at 1 T it is nearly an or- 
der of magnitude higher and starts to increase strongly 
with decreasing temperature, see left inset in Fig. [I] The 
nature of the gap opening at the CNP is elucidated fur- 
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FIG. 1. Resistance as function of the concentration by sweep- 
ing the backgate from -60 V to 60 V for T (black ■ ) and 
1 T (red •). A constant contact resistance has been re- 
moved, top-left inset: Temperature dependence of the 
resistance at the CNP for T and 1 T; top-right inset: 
SEM-picture of our suspended device. 



ther in Fig. [2^, where we show the resistance as a func- 
tion of carrier concentration n for several magnetic fields. 
The diverging resistance at the CNP appears at similar 
magnetic field as the plateaus at filling factors -4 and 
4; i. e. the eight-fold degeneracy of the zero-energy Lan- 
dau level breaks directly into two four-fold degenerated 
Landau levels, as already predicted theoretically^ and 
proven experimentally.^ At low fields B < 0.1 T we ob- 
serve a small decrease of the resistance maximum at the 
CNP (not shown in the figure). This small decrease in 
resistance can be explained by the presence of local in- 
homogeneities which give a small splitting between the 
valley-polarized energies; the cross-over of these energy- 
states at finite magnetic field results in a resistance min- 
imum. When the magnetic field is above B > 0.1 T we 
observe a rapid increase of the resistance-maximum at 
the CNP, shown in Fig. [2|d. We can interpret this rapid 
increase as a result of the spin-splitting of the two en- 
ergy levels at zero energy or by disorder, e.g. unevenly 
charged top and bottom layer. The last scenario would 
lead to a strong temperature dependence at zero field and 
ultimately for big disorder to an insulating state at zero 
field, as discussed in Ref. [28] The absence of a tempera- 
ture influence at T and the CNP centered at very low 
gate- voltage points to a non-disordered bilayer, therefore 
we interpret the rapid increase by a result of spin split- 
ting. 

The absence of an energy-level at E = in the in- 
set of Fig. ^jp results in a diverging resistance at the 
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CNP. The resistance Rcnp at the CNP follows a classi- 
cal Ahrrenius- activation behaviour R xx oc exp (A/fc^T), 
in where A is a scale for the size of the gap. The re- 
sistance increase scales best with In (R) oc B /T, from 
which we obtain a gap A = 0.34 meV/T x B. This 
gap is about a factor 3 times larger than the Zeeman- 
splitting gn B B, which can be explained by the dominat- 
ing exchange energyP^ Equation ([lj describes the total 
spin energy As, determined by the sum of the single 
electron Zeeman energy gfisB and the exchange energy 
E ex ' ( n -f ~ n i)- Here — is the normalized difference 
between spin-up and spin-down occupation. 



A s = 9VbB + E ex ■ (n t - nj.) 



(1) 



At low fields the two energy levels are still overlapping 
and the system is not fully spin polarized, (n^ — n^) < 1. 
Assuming Gaussian shaped Landau levels we can approx- 
imate (n^ — n-i) = „^r- with leads with help of equa- 
tion (jlj) to the gap Ag = 1 f^ B ^ r - The observed spin- 
enhancement by a factor 3 corresponds to a typical level 
width r = 2 meV and exchange energy E ex = 1.3 meV 
at B = 1 T corresponding to a value of about 2 % of the 
Coulomb energy E c — e 2 /e r lB — 56 meV, where l B is 
the magnetic length. 

The behavior at high magnetic field is experimentally 
more complicated to access, because the measured re- 
sistance rapidly exceeds several Mf2s and a quantitative 
analysis becomes difficult. However away from the CNP 
the measured resistances stay low enough to guarantee a 
reliable interpretation up to the highest magnetic fields. 
This situation is illustrated in the inset of Fig. [2Jd where 
we sketch the quantized density of states in the lowest 
Landau level around the CNP with a gap 2A opening 
at E — 0. When the Fermi energy is located at a finite 
energy E < A (i.e. still inside the localized parts of the 
DOS), conduction will occur by thermal excitation to the 
conductivity edges E c \ and E C 2 of the extended states. 
The resistance R(E) at this energy will then be given by 



A(E) + E A{E)-E 

R(E) oc e k » T +e k » T 



1 A(g) 

-e k <> T 
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E 
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For relatively small energies E « k^T the cosh-term 
can be approximated by a first order Taylor expansion 
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7(.E). For small E we 



can interpret equation (2 1 as Rcnp oc \e k >> T j(E). At 
the CNP, E — 0, this approaches a trivial Ahrrenius- 
behavior, while for non-zero fixed energy ^y(E) is an 
energy dependent renormalization factor which for 
E « kbT is independent of T. 

We analyze the resistance at concentrations 
n = ±3 x 10 14 m~ 2 (dots marked in Fig. [2^, and 
multiply this data with a fixed constant to make 
an overlap with the low field data). All datapoints 
R > 1 Mf2 in Fig. are verified by this method and 
therefore reliable up to the highest field. This proper 



scaling for both low and high resistances also excludes a 
strong effect of the local heating due the finite excitation 
voltage we apply over the sample. 

From Fig. [2]d we see that the scaling of the resistance 
at high magnetic fields is remarkably different from the 
linear field increase at low fields. This observation is 
again visualized in Fig. Ep where we show that the 
calculated gap strongly bends and the slope strongly 
reduces. In this regime the gatesweeps are packed more 
densly for increasing magnetic field and the energy E 
gets comparable to the thermal activation k^T thus 
we are no longer able to calculate the gap-size with a 
simple Arrhenius-behaviour. Experimental limitations 
of our suspended samples do not allow us to access much 
higher temperatures, therefore we can only speculate 
here about further gap-study. 

At high enough fields we expect to fulfill the criteria of 
fully spin-polarized system, (n-f — n\) — > 1. The sudden 
strong change of the gap-size suggests that our system 
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FIG. 2. (a) Resistive behaviour of the sample near the CNP at 
T—1.3 K. The dots mark the positions n = 3 x 10 14 , and — 
3 x 10 14 m~ 2 where the gap opening has been analyzed (see 
text for more details), (b) Resistance Rcnp as function of 
B/T at T — 4.2 K and 1.3 K; inset: Qualitative picture of 
the DOS near the CNP. The conduction at energy E is di- 
rectly related to the thermal excitation of electrons to the 
conduction edges E c i and E C 2- (c) Calculated gap A as func- 
tion field B for T = 1.3 K; the dashed line represents the 
theoretical single electron Zeeman-energy q^lbB. 
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FIG. 3. (a) Gatesweeps R(Va) at constant magnetic fields 
B = 1 T, 5 T, 12 T, 17.5 T and 30 T for T = 4.2 K; the 
curves are shifted up for clarity, (b) Derivative dR/dVG for 
the curves in (a), (c) Position of the minima for v = ±12, 
±8, ±4, ±3, ±2 and ±1 as function of the magnetic field. 



indeed gets fully spin-polarized, in literature also known 
as the cross-over to a quantum Hall ferromagnetic 
state. Further increase of the magnetic held leads 
hypothetically to a dominating spin-splitting gfisB, 
because E ex oc \AB. Further specific research in titled 
magnetic helds is necessary to decouple the influence of 
the single electron Zeeman energy and exchange energy. 

After detailed study at low concentrations we have a 
closer look at the manifestation of the QHE at higher con- 
centration. In Fig. [3^ we show the corrected two-point 
resistance as a function Vq at 4.2 K for B = 1 T, 5 T, 
12 T, 17.5 T and 30 T. Apart from the distinct v = ±4 
plateaus which are already well pronounced at 1 T, ad- 
ditional plateaus at v = ±3, ±2, and ±1 start to appear 

in higher fields. In Fig. J^Jd we show the derivative 

of the resistance curves, where we can already already 
recognize distinct maxima and minima for lower fields. 

In Fig. [3J; we follow the position of the minima with 
increasing magnetic field. We see that the maximum 
applied gate-voltage Va = 60 V limits the observation 
of filling factors v — ±4 up to 9 T, while v = ±2 re- 



FIG. 4. (a) Tia^-oscillations after removing linear background 
from n|f for B = 3 T, 5 T, 12 T and 20 T at 4.2 K. (b) Din- 
gle plot of v — 8 and v = 4: amplitude A of the oscillations as 
function of the inverse field 1/B. (c) Dingle plot for v = 3, 2 
and 1. (d) Schematic plot of the appearance filling factors 
with increasing magnetic field. The dashed lines show the 
center of the Landau level, while the grey shaded area is the 
Landau level broadening determined by the Dingle tempera- 
ture To- 



mains observable till fields of 15 - 20 T and filling factor 
v = ±1 is still observable at the highest applied field, 
30 T. From Fig. [3]: we observe that the position of the 
minima strongly deviates from the linear relationship be- 
tween the induced charge carrier concentration n and 
the applied magnetic field B, n = v^r-. The equidis- 
tance of the minima for fixed magnetic fields excludes 
a capacitance-change due the bending of the membrane; 
which can be expected by the particular big difference 
between the electric field induced bending (10 - 20 nm) 30 
and the vacuum gap over which graphene is suspended 
(~ 1.5 /im). In the appendix we discuss in more detail 
how to extract the exact relation between the applied 
field B and the induced charge concentration n from this 
data. 

Experiments on suspended graphene-samples are 
mainly performed in two-probe configuration. Ex- 
perimental limitations of the annealing-procedure do 
not allow us to obtain very homogenous samples in 
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four-probe configuration. Therefore more effort has to 
be done to do a proper analysis on both the magnetore- 
sistance and Hall-resistance. In Fig. [4] the appearance 
of the different filling factors are further elucidated; 
in particular at positive gate-voltages, the influence of 
contact-resistance is here experimentally the smallest. 
As shown in Fig. |3jb the derivative of our data, jp^, 
shows already at very low field a very clear appearance 
of filling factors v = 1, 2, 3 and 4. A small change 
in the slope of R xy causes a very distinct minimum in 
the derivative. Theoretically we can use a model that 
directly describes the magnetoresistance R xx in terms 
of the Hall resistance R xy W^ i- e - Rxx oc n d ^ y . In 
Fig. |4^l we study the appearance of the filling factors 
by plotting the obtained magnetoresistances R xx — Rq. 
Here we removed from all data the linear background 
Ro of the 12 T data shown in the inset of Fig. [4^, 
and centered all curves around the x-axis. We used 
the obtained leverage factor a(B) to determine the 
exact concentration n. Already at 3 T we observe 
the appearance of clear oscillations around v = 2 and 
v = 3 followed by the appearance of v = 1 at 5 T. 
The amplitude A of the oscillation is defined by the 
difference between the minimum and the first maximum. 
A single oscillation can be best analyzed by applying 
Lifshitz-Kosevic equation^ A ■ cos (f(B)), here A is the 
amplitude and f(B) a field-dependent function that 
determines the frequency and phase. The amplitude 
A is finite due to the Landau-level broadening, and 
is damped by the Dingle-factor exp (— /3 ■ T£>m c /B) in 
where To is the Dingle-temperature, m c the cyclotron 
mass in units of electron mass m e , B the magnetic 
field and (3 = 14.694 T/K. In Fig. [Ija the amplitudes 
A for v = 4 and v = 8 are plotted as function of 1/B, 
which affects in a linear decrease with slope j3 ■ Tj^m c . 
If we assume the cyclotron mass in bilayer graphene 
to be m c fti 0.033 • m e (corresponding to 71 = 0.4 eV, 
see Ref. [33J and references in there) we obtain the 
Dingle temperatures Td in the table below. We repeat 
the same procedure for filling factors 1, 2 and 3 in Fig.|4j;. 



V 
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Td (K) 


2.4 ±0.4 


1.4 ±0.4 


29.2 ±4 


9.2 ± 1.2 


58 ±8 



Compared to v = 4, fully quantized at B = 1 T, the 
Dingle temperatures Td for the degenerate filling factors 
are one order of magnitude larger, which means fields 
B > 10 T are required to observe full quantization; in 
particular filling factor v = 1 becomes quantized at fields 
of B > 30 T. 

In Fig. [4] we illustrate qualitatively the appearance of 
the different Landau-levels for increasing magnetic field. 
The corresponding grey shaded areas describe the Lan- 
dau level broadening T, directly determined by the Din- 
gle temperature To; higher Dingle temperatures corre- 
spond broader Landau levels. While the position of the 
energy moves linearly with increasing field, the Landau 
level broadening T is proportional to the square root of 



the applied field r oc \[B. With increasing field the over- 
lap between the shaded areas decreases, and the plateau 
starts to appear. As we can see from Fig. |4ji Landau 
levels around v = 2 and v = 3 do indeed not overlap 
anymore for similar magnetic field, however the shaded 
areas for v — 1 overlap till higher fields. Finally the 
overlapping of filling factors v = ±1 disappears at simi- 
lar magnetic field as the resistance at the CNP starts to 
bend strongly, supporting the idea of a cross-over to a 
fully spin polarized state at v = 0. 

After the appearance of the non-degenerated filling fac- 
tors v = 4 and 8 a gap at v = forms, followed by filling 
factors v = 2 and at high fields v = 1 and 3. This se- 
quence agrees with the proposed model of the formation 
of a quantum Hall ferromagnet in the lowest Landau level 
and the observed behaviour at the CNP. 



IV. CONCLUSION 

In conclusion we have performed experiments on a 
suspended BLG sample which shows us a field induced 
gap at the CNP for fields B > 1 T. The gap at v = 
opens simultaneously with the formation of filing factors 
v = ±4, which implicates the eight-fold degenerated low- 
est Landau breaks directly in two fourfold-degenerated 
spin-polarized subbands. At high magnetic fields we 
observe a smooth transition to a much smaller gap, 
this is consistent with the picture of the formation 
of a spin-polarized quantum Hall ferromagnetic state. 
Following the breaking of the lowest Landau level we 
observe a breaking of v = in v = ±2 and finally in 
v = ±1 and v = ±3, in agreement with the theoretically 
proposed model of a quantum Hall ferromagnet. 
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V. APPENDIX 

Fig. [5] shows the relation between a and the applied 
magnetic field. The capacitance of the sample increases 
from the geometrical value 0.5 xlO 14 m~ 2 V _1 up to 
1.8 xlO 14 tot 2 V _1 at 9 T and saturates at this value 
for the highest fields. Thi s eff ect is also observed implic- 
itly in recent publications^^ on high quality suspended 
bilayer-devices, but not mentioned by authors in the text. 
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As discussed in ref. the increase in capacitance of 
the system under applied magnetic field could be under- 
stood from the deviation from the flat-plate capacitor 
model. At the point when the width of the graphene 
flake is smaller or comparable to the distance to the back 
gate the flat-plate capacitor model can be no longer ap- 
plied. The charge carrier distribution in graphene be- 
comes non homogenous and increases at the edges. Since 
the classical cyclotron radius of charge carrier depends 
inversely on magnetic field, the increase of B will cause 
edge channels in quantum Hall regime propagate closer 
to the edge, where the density can be few times higher 
than in bulk graphene. This would result an increase 
in capacitance extracted from QHE plateaus. The cy- 
clotron radius is expected to be dependent on the charge 
carrier density as well. The exact calculations for differ- 
ent device's geometries with charge carrier distribution in 
graphene, compared to the experiments, are the subject 
of another paperf^S 
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FIG. 5. Relation between the induced charge carrier concen- 
tration per applied magnetic field 01(B) = n/B and the ap- 
plied magnetic field B for B=0 (•) and B (o). The plotted 
line shows the interpolated a(B) from which we determined 
a reliable value of the concentration n. 
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